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Abstract. - We studied interplay between kinetic roughening and phase ordering in 1+1 
dimensional single-step solid-on-solid growth model with two kinds of particles and Ising-likc 
interaction. Evolution of both geometrical and compositional properties was investigated by 
Monte Carlo simulations for various strengths of coupling. We found that the initial growth 
is strongly affected by interaction between species, scaling exponents are enhanced and the 
ordering on the surface is observed. However, after certain time, ordering along the surface 
stops and the scaling exponents cross over to exponents of the Kardar-Parisi-Zhang universality 
class. For sufficiently strong strength of coupling, ordering in vertical direction is present and 
leads to columnar structure persisting for a long time. 



Recently there has been considerable interest in growth induced surface rough- 
ening called kinetic roughening and in phase ordering independently. 
In kinetic roughening one is interested in the evolution of roughness during 
a nonequilibrium growth process. Phase ordering deals with the approach to 
equilibrium of a system quenched from a homogeneous high temperature phase 
into a two-component region. In both fields concept of scaling allowed to classify 
physical processes into universality classes. However, not very much is known 
so far about scaling exponents in systems where both processes, roughening as 
well as ordering, are relevant. 

Growth of many-component systems is in fact quite common situation. Many 
real materials are composed of two or more components but often the dynamics of 
one component is dominant and one can use a single- component growth model. 
For example, one can consider only kinetics of Ga atoms in studies of GaAs 
growth ]3|]. We are interested here in a more complex case when dynamics of 
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both components is important. It is a problem of practical interest, for example, 
microscopic understanding of growth of alloys is desirable (|, ||, ||. One can study 
different aspects: kinetics, morphology, scaling behaviour etc. We are not going 
to model growth of any specific material but we shall concentrate on the scaling 
behaviour. The problem of growth in a system with two or more components 
is interesting from pure statistical-mechanical point of view, because growth 
process may belong to a new universality class [0, |^, [J ; such system might also 
exhibit a nonequilibrium phase transition between the low and high temperature 
region. 

Let us consider a surface in a ci-dimensional space given by a single valued 
function h(r, t) of a cf-dimensional (d=d'+l) substrate coordinate r. The surface 

roughness is described by the surface width w(t,L) = (y h? — h?), where t is 
the time, L is the linear system size and the bar denotes a spatial average, 
(...) a statistical average. The surface roughness often obeys dynamical scaling 
law w(t, L)(xL^ f [t / L z ) where the scaling function f(x) has properties: f(x)= 
const., 2^>1 and f(x)oax 13 , x<^l {(3—Q/z). The exponents ( and z (or ( and 
f3) characterize scaling behaviour of the surface width in a particular model and 
determine its universality class [|]]. This universal behaviour has been observed in 
a wide variety of growth models and there has been considerable effort in finding 
different possible universality classes. Many of the growth models studied so far 
(for example ballistic deposition, Eden model, restricted solid-on-solid model, 
etc.) belong to the Kardar-Parisi-Zhang (KPZ) universality class [ID . 



There is a large menagerie of single- component growth models which can 
be potentially generalized to heterogeneous case. Moreover, there are different 
possible ways of generalization. Nevertheless, little is known so far about kinetic 
roughening in two-component growth models. This problem was probably first 
considered by Ausloos et al. J7j. They introduced a generalization of the Eden 
model coined as magnetic Eden model (MEM), which contains two types of 
particles with probabilities of growth given by Ising-like interaction. Ausloos 
et al. found a variety of morphologies, and they also measured the perimeter 
growth exponent 0. They suggested that the model does not belong to the KPZ 
universality class. Recently Wang and Cerdeira studied kinetic roughening 
in two 1 + 1 dimensional two-component growth models with varying probability 
of deposition of a given particle type. They did not found, however, a new 
universal behaviour. Although the phase ordering was apparently present it was 
not studied in these works. The characteristic length in phase ordering is a 
domain size D. It increases with time according to a power law, D oc i*, i) 
being one of the exponents specifying a universality class. Phase ordering is 
usually a bulk process, but growth induced ordering may be present just on the 
surface. Then evolution of the domain size on the surface is of interest. The 
scaling in this case has been recently studied by Saito and Muller-Krumbhaar 
H in another generalization of the Eden model (different from MEM). They 
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obtained ip — 2/3 for growth of domain size, but they did not investigate kinetic 
roughening. 

In this Letter we study interplay between kinetic roughening and phase order- 
ing. We introduce a new growth model with two types of particles suitable for 
the study of scaling and present results of simulations in 1+1 dimensions. Our 
model is based on the single-step solid-on-solid (SOS) model, a discrete model 
with constraint that the difference of heights between two neighbouring sites 
is restricted to +1 or —1 only. Advantage of this choice is that in contrast to 
cluster geometry, kinetic roughening can be more easily studied and that the 
single-step constraint allows simple interpretation of the deposition process in 
a binary system. We consider two types of particles, and denote the type of 
a particle by a variable a which assumes values +1 or —1 (not to be confused 
with the step size). It is convenient to use an analogy with magnetic systems 
and to consider a as a spin variable. We shall use this terminology in the 
following, but it should be warned that in the context of crystal growth this may 
be misleading because magnetic interactions of atoms are rather weak and do 
not usually control the growth process. One should rather think in terms of two 
types of particles with different bonding energies in this case. Therefore we shall 
call our model two- component single-step (TCSS) model. 

We describe our growth model for simplicity in 1+1 dimensions but it can 
be straightforwardly generalized to any dimension. Several realization of the 
single-step geometry differing by the number of nearest neighbours for a new 
particle are possible. Here we consider a variant with three nearest neighbours 
which can be represented as stacking of rectangular blocks with the height twice 
the width. During growth, particles are only added, there is neither diffusion 
nor evaporation. Once a position and a type of a particle are selected, they 
are fixed forever. Due to the single step constraint, particles can be added 
only at sites with a local heigh minimum, called growth sites. The probability 
of adding a particle with a spin a to a growth site % depends only on its 
local neighbourhood in a way analogous to rules used in the magnetic Eden 
model 0. It is proportional to exp(—AE(i, a) / ksT) , where ks is Boltzmann's 
constant, T, and AE(i, a) denote thermodynamic temperature and change of 
energy associated with deposition of a new particle, respectively. The energy 
AE(i, a) is given by Ising-like interaction of a new particle with particles on 
the surface within nearest neighbours of a growth site (which are three in the 
selected realization - left, bottom and right). When we denote the spin of a 
particle on the top of a column of spins at site i (surface spin) by o~(i) then 
AE(i, a) = — J a [a(i — 1) + a(i) + a(i + 1)] — Ha. Here J is a coupling strength 
and H is an external field. In the following we use for convenience dimensionless 
constants K = J/ksT and h = H/ksT. 

We studied both geometrical and compositional properties. Geometry is de- 
scribed by the surface width w(t, L) or by height-height correlation function 
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G(r,t) = j^J2i=i([h{i + r,t) —h(i,t)] 2 ). Evolution of geometry is affected by 
composition of the surface and vice versa. Therefore we introduce quantities 
characterizing composition on the surface. We concentrate only on ordering on 
the surface, we do not study bulk properties here. We consider the average 
domain size along the surface D(t), the spin correlation function of surface spins 
S(r,t) = |EiiWj + f,<WM)>, and magnetization M(t) = x Ef=i(o~(M))- 
We performed simulations for various coupling strengths K in both ferromagnetic 
(K > 0) and antiferromagnetic (K < 0) regimes, mostly for zero external field h. 
System sizes ranged from L = 250 to L — 80000, time was up to 3 x 10 5 monolay- 
ers (ML). Presented results are averages over ten or more independent runs. We 
start from the flat surface as usual, but in two-component models the evolution 
strongly depends on initial composition of the substrate. We considered various 
possibilities (ferromagnetic, antiferromagnetic, random composition) and finally 
we have used growth on a neutral substrate, i.e. a substrate composed of particles 
with zero spin, and we let the system to evolve spontaneously in order to avoid 
initial transient effects. 

Fig. 1 shows examples of time evolution of morphologies obtained for selected 
couplings. We can see that with the increasing coupling the surface becomes 
more and more rough (facetted) and at the same time larger and larger domains 
of the same type of surface particles are formed. Notice also that for the larger 
coupling there is correlation between domain walls and changes of the local slope, 
and a columnar structure is observed. The time dependence of the surface width 
and the average domain size for L = 10000 are shown in Fig. 2 and Fig. 3, 
respectively. When the coupling is weak evolution of the roughness is almost 
the same as in the ordinary single-step model (f3 = (3^ KPZ ^ = 1/3). For a larger 
coupling it can be seen that the average surface width at a given time is an 
increasing function of the coupling and that ordering leads to an enhancement 
of the exponent /5 e ff ■ However, after a time t cross , that increases with the coupling 
strength, the exponent (3 eS crosses over back to f3^ KPZ ^ = 1/3. To be sure that 
this unexpected crossover is not a finite size effect we performed for K — 1 
the calculation for significantly larger system size L = 80000 and observed the 
same crossover (see inset in Fig. 2). When the coupling is even stronger the 
enhanced exponent (3 e g is observed during the whole simulation; for K = 2 we 
have (3 cS = 0.52 ± 0.02 for times up to 3 x 10 5 ML. 

The crossover can be understood from behaviour of the average domain size. 
It increases initially according to a power law D oc t^ cS , ip e ff being an increasing 
function of K; ip e g = 0.33 for K = 1.1. However, after a certain time t^t , D 
saturates to D sat (K). We have checked that when the saturation is observed 
it is not a finite size effect (see inset in Fig. 3), rather an intrinsic property 
of the model. For large K the domain size increases during whole simulation, 
after some transient effect we observe ifj e g = 0.44 for K = 2. The calculation 
of the spin-spin correlation function, and the correlation length £ s derived from 



M. kotrla et ai: INTERPLAY between kinetic etc. 



5 



it, leads to similar results: during ordering the correlation length increases as a 
power law £ s (t,K) oc t KcB ; K e g « 1/2 for K = 2. Then after a time increasing 
with the coupling ordering stops, and £ s saturates to £| at oc e yK ; 7 = 3.25 ±0.08. 
Both effects, the crossover in time dependence of roughness and the saturation 
of the domain size (correlation length), are apparently related, t cross ~ (see 
inset in Fig. 2). 

To complete the picture we need to know the second exponent for kinetic 
roughening, the roughening exponent £. We measured it from the spatial depen- 
dence of the height-height correlation function G (r,i) at £g™ t (not shown here). 
The exponent has a value close to Q KPZ > = | for a weak coupling, but for a 
larger coupling we have found that there is a crossover behaviour provided the 
system is sufficiently large. We observed that the exponent £ e ff is increasing with 
the coupling (£ e ff = 0.7 for K = 0.7) on distances smaller than the correlation 
length £ s , and c^ e ff crosses over to ^ KPZ > = ~ on a larger distances. If the coupling 
is strong we again do not see the crossover but only larger exponent, Cofr ~ 1 
for K = 2, because the simulated system sizes are not large enough to get into 
the regime r > £ s . Exponent ( = 1 is the same as in some models with surface 



diffusion in which often anomalous scaling is observed ||11|| . However, this cannot 
be the case here since the step size is restricted. 

Our results show that for sufficiently strong coupling the TCSS model exhibits 
an intermediate growth regime with new scaling exponents which we estimate 
as/5 = 1/2, ( = 1 (z = 2) and ip = 1/2. However, these exponents are 
only effective and there is a crossover to the KPZ exponents for the medium 
coupling strengths. We estimated the dependence of t^t (~ Across) on K. We 
have found that the saturated domain size as a function of K can be well fitted 
in the form D SSut (K) = 1 + |(e 35X + e 0JK ) (see inset in Fig. 3). Comparing 

D BS t (tgS )^ with the fit we obtain that t^t is a rapidly increasing function of 
K, t^(K) oc e 7K , e.g. t^(K — 2) w 10 6 . Due to the progressively increasing 
time and system sizes needed for simulation we cannot exclude that there is a 
phase transition to a new true asymptotic phase at some value K c . However, 
we expect that the crossover in the roughness is present for any value of K, 
but that it is hard to see it for strong coupling because t cross is larger than any 
possible simulation time. When K = 00 then the model becomes trivial: type 
of a new particle added to a growth site is dictated by the majority of particle 
types in the neighbourhood and evolution of composition is fully determined by 
the initial configuration. So far all results were for zero external field. Non-zero 
field leads to nonzero surface (as well as bulk) magnetization, or in the context 
of alloy growth to changing stoichiometry. We can still define the exponent 
ijj for growth of the dominant domain size as well as the exponents for kinetic 
roughening. Effect of this symmetry breaking on the values of exponents remains 
to be studied. 
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One would like to know if the crossover to the KPZ exponents is generic 
behaviour or if facetted phase can be present also in the asymptotic regime. We 
have also tested a different variant of two-component growth, namely the single- 
step model which corresponds to stacking of squares rotated by 45 degrees in 
which a new particle interacts with only two instead of three nearest neighbours. 
We observed the crossover as well, stronger coupling was needed to see enhanced 
exponents and for the same coupling the crossover time was smaller than in 
the variant with three nearest neighbours. Is is of interest to examine kinetic 
roughening and phase ordering in a different two-component growth model, in 
particular to reexamine the MEM model in the strip geometry since in the 
cluster geometry used by Ausloos et al. finite size effect are expected to be 
stronger. Situation in higher dimensions can be investigated by straightforward 
generalization. Finally, there is question of continuous description of the two- 
component growth. 

In conclusion we have generalized the single-step SOS model to two-component 
growth model. It can be considered as a model for growth of binary alloys 
from a fluid containing two types of particles, or for growth of a colony of two 
kinds of bacteria etc. The introduction of two kinds of interacting particles 
leads to new phenomena, like an increase of the roughness with the increasing 
coupling strength, larger effective scaling exponents and ordering on the surface 
with corresponding pattern formation during growth. However, after a certain 
time new behaviour stops and ordinary dynamic behaviour is restored. Hence, 
growth of a crystal can be divided into two stages. In the first, ordering is 
essential and influences the evolution of geometric properties. In the second 
stage, after saturation of domain size, the evolution of geometric properties is 
similar to that of the ordinary single-step model. Crossover time between two 
regimes is independent of system size and is a rapidly increasing function of 
coupling so that practically only the first regime may be observed. Patterns of 
the grown material in saturation regime are similar to lamellar structure observed 
in eutectic growth J|. At the moment, it is not clear why ordering stops and 
the characteristic length (given by the saturated domain size) is selected. 
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Fig. 1. - Examples of evolution of surface profiles for several values of coupling strength, 
K = 0.3, 0.7, 1.1, 2.0 and zero external field. Surface profiles at various times increasing 
as powers are shown by white lines. Only part of the substrate close to the surface is shown at 
given time, black and grey correspond to different types of particles. System size is L — 250. 



Fig. 2. - Surface width w vs. time t for several values of coupling strength, K = 0.0, 0.3, 0.7, 
1.1, 2.0 and zero external field, L = 10000. Inset: Comparison of the time dependence of the 
surface width and the domain size for K = 1.0 and L = 80000. 



Fig. 3. - Time evolution of the surface domain size for coupling strengths K — 0.0, 0.3, 0.7, 
1.1, 2.0 and zero external field, L — 10000. Inset: Saturated surface domain size as function 
of coupling for different system sizes L — 250, L = 1000 and L = 10000. The solid line is the 
fit (see text). 
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